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(Dated: February 7, 2008) 

We present a search for anomalous production of heavy-flavor quark jets in association with a 
W boson at the Fermilab Tevatron pp Collider. This search is conducted through an examination 
of the exclusive jet spectrum of VF+jets final states in which the heavy-flavor quark content has 
been enhanced by requiring at least one tagged jet in an event. Jets are tagged by the combined 
use of two algorithms, one based on semileptonic decays of b/c hadrons, and the other on their 
lifetimes. We compare data in e + jets (164 pb _1 ) and fi + jets (145 pb _1 ) channels, collected with 
the D0 detector at *Js = 1.96 TeV, to expectations from the standard model, and set upper limits 
on anomalous production of such events. 



PACS numbers: 13.85.Qk, 13.85.Ni 



The heavy-flavor (HF) content of jets produced in as- a non-SM source of physics. The CDF collaboration rc- 
sociation with a W boson in pp collisions provides a test cently reported just such an excess in the exclusive W+ 
of the standard model (SM), and an excess would suggest HF-jet spectrum in which one jet was tagged using both 
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secondary-vertex (SVT) and soft-lepton (SLT) tagging 
algorithms . To check for the presence of this anomaly 
in our data, we also select jets tagged with both algo- 
rithms. In addition, we use two benchmark SM processes 
as models for new physics and derive upper limits on such 
processes. 

At the Tevatron, the primary SM contributions to a W 
boson associated with HF quarks in the final state are ex- 
pected to be from tt, Wbb/cc (where the bb or cc pairs 
arise from gluon splitting), and Wc final states, with ad- 
ditional contributions arising from single top quark or 
WZ (with Z — > bb/cc) production. The production of W 
bosons accompanied by light quarks or gluons (referred to 
as VF+jets in this Letter) contributes to the background 
when the light-quark or gluon jets are misidentified as jets 
from HF quarks. Since W bosons are identified through 
their W — > ev and W — ► [iv decays, background can arise 
from Zbb, ZZ (with one Z — ► bb/cc), and Z+jets produc- 
tion when one of the leptons from the Z — >■ £ + £~ decay 
is not observed in the detector. The main instrumen- 
tal background arises from multijet processes in which 
a jet is misidentified as a lepton, and an imbalance in 
transverse momentum {I/t) is generated through a mis- 
measurement of the jets or a lepton. To be selected, 
these kinds of events must also contain tagged HF jets or 
misidentified non-HF jets. 

The data were collected with the D0 detector 2] dur- 
ing Run II of the Fermilab Tevatron Collider in pp col- 
lisions at y/s = 1.96 TeV. The components used in 
this analysis include the central tracker, calorimeter, 
and muon detectors. The central tracker consists of 
a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT), both located within a 2 T supercon- 
ducting solenoidal magnet. The uranium/liquid-argon 
calorimeter consists of three sections, each housed in a 
separate cryostat The central calorimeter (CC) cov- 
ers pseudorapidity \rj\ < 1.1, while the two end calorime- 
ters (EC) extend the coverage to \r)\ w 4.0. The muon 
system is located outside the calorimeters, and consists 
of a layer of tracking detectors and scintillation trigger 
counters inside 1.8 T iron toroids, followed by two more 
similar layers outside the toroids. 

The W — + eu and W — > \iv decay candidates are 
selected initially by triggering on electrons and muons. 
The average trigger efficiency for electrons with trans- 
verse momentum px > 20 GeV/c and \r)\ < 1.1 is 
(97.0±0.3)%. The average trigger efficiency for muons 
with p T > 20 GeV/c and \r)\ < 1.6 is (62.1±3.4)%. The 
integrated luminosity is 164 ± 11 pb _1 for the electron 
sample and 145 ± 9 pb _1 for the muon sample. 

Candidate events for W — ► eu decays are selected 
by requiring exactly one isolated electron with px > 
20 GeV/c and < 1.1, defined relative to the geometri- 
cal center of the detector. Lepton isolation requires a sep- 
aration in rj and azimuth (</)) of TZ — v/ (Ary) 2 + (A0) 2 > 
0.5 from all jets in the event. Electrons are defined us- 



ing a cone algorithm, and by the energies deposited in 
calorimeter towers within a radius of TZ = 0.2 of the elec- 
tron axis, with at least 90% required to be within the 
electromagnetic portion of the calorimeter, and by the 
total energy in a cone of TZ = 0.4 centered on the same 
axis, which must not exceed by more than 15% the recon- 
structed electron's energy. In addition, the longitudinal 
and transverse shower shape must be compatible with 
that expected from an electron. 

Candidate events for W — > [iv decays must contain 
exactly one isolated muon with pr > 20 GeV/c and \rj\ < 
1.6, also defined relative to the geometrical center of the 
detector. Muons are required to satisfy two additional 
isolation criteria: the transverse energy deposited in the 
calorimeter in the annular region of 0.1 < TZ < 0.4 around 
the muon's path must be smaller than 2.5 GeV; and the 
vector sum of the px values of all tracks within TZ = 0.5 
of the muon's trajectory must be less than 2.5 GeV/c 
(excluding the track matched to the muon's trajectory). 

Lepton identification is refined by requiring the trajec- 
tory of a track reconstructed in the SMT and CFT to 
match either the position of the electron energy cluster 
in the calorimeter or the position of hits in the muon 
detector. To complete the selection, all events are also 
required to have I/t > 20 GeV, and the azimuthal an- 
gle between the lepton and the direction of the Jfx must 
be greater than 7r/8. To eliminate poorly reconstructed 
events, the primary vertex (PV) of the event must con- 
tain at least three tracks, and its z-position (along the 
beam) has to be closer than 60 cm from the center of 
the detector. Finally, to reject multijet background, we 
require a reconstructed transverse mass consistent with 
that of the W boson, 40 < M Wt < 120 GeV/c 2 . In cal- 
culating Mw T , we assume that the IpT corresponds to 
the transverse energy of the neutrino. 

Upon selection of VF-boson candidates, we evaluate the 
HF-quark content of each event. Jets are defined using 
an iterative seed-based cone algorithm (including mid- 
points), clustering calorimeter energy within TZ = 0.5. 
This is subsequently corrected for jet energy scale, based 
on momentum balance in photon+jet events . We con- 
sider only jets with Et > 25 GeV/c and |?7| < 2.5. These 
jets are then evaluated using two HF-tagging algorithms, 
as described below. 

The soft-lepton tagging (SLT) algorithm is based on 
low-px muons arising from semileptonic decays of HF 
quarks (via virtual W bosons) that are produced near 
a jet in (j], 4>) space. Only muons with pt > 4 GeV/c 
and 1 17 1 < 2.0 are considered. To reject Z — > /i/i back- 
ground, we require pt < 15 GeV/c for the muon. Jets 
with a muon within TZ = 0.5 of the jet axis are deemed 
tagged. Typical SLT efficiencies for &-quark jets are ap- 
proximately 11%, and 0.4% for light-quark jets. The 
additional muon present in SLT events causes an in- 
crease in the average single-muon trigger efficiency from 
(62.1 ± 3.4)% to (68.4 ± 3.5)%. 
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Secondary- vertex tagging (SVT) is used to identify dis- 
placed decay vertices of long-lived particles. To form 
secondary vertices (SV), charged tracks are selected on 
the basis of the significance of their distance-of-closest- 
approach (dca) to the PV. Tracks are first grouped in 
1Z = 0.5 cones around a seed track with pr > 1 GeV/c 
and dca/ddca > 3.5, where Odea is the uncertainty on the 
track's dca. Proto-vertices are formed by adding tracks 
to the initial grouping, provided their contribution to the 
X 2 of the vertex fit is small. Secondary vertices are se- 
lected by requiring the transverse distance from the SV 
to the beam direction, L xy , to be less than 2.6 cm, and 
the decay-length significance, L * v to be greater than 7, 
where cr^ xH is the estimated uncertainty on L xy calcu- 
lated from the error matrices of the tracks in the vertex. 
Jets are considered tagged by this algorithm when a SV 
lies within 1Z = 0.5 of the original jet axis. This SVT 
algorithm exhibits a typical tagging rate for 6-quark jets 
of ~ 32%, and 0.25% for light-quark jets. 

To predict SM rates, Monte Carlo (MC) events are 
generated for the processes mentioned above, with the 
exception of multijet production, which is estimated from 
data as described below. W/Z+jets (both HF and light- 
quark jets), tt, and diboson processes are simulated 
with ALPGEN 5]. Single top quark processes are sim- 
ulated using COMPHEP p . All events are generated with 
mtop = 175 GeV/c 2 . Hadronization and showering of 
these events is based on pythia [tJ. The exceptions are 
W/Z + b processes, where Zb is simulated using PYTHIA, 
and the contribution from Wb is estimated from the pa- 
rameterized mcfm MC Q and used to calculate a cross 
section relative to Wbb production assuming the jet-pr 
spectrum from pythia for inclusive W boson produc- 
tion. All MC events are generated at y/s — 1.96 TeV, 
using CTEQ5L parton distribution functions and a de- 
tailed detector simulation based on GEANT Ifj. To simu- 
late the effect of multiple interactions in beam crossings, 
a Poisson-distributed minimum-bias event overlay, with 
an average of 0.8 events, is included for all events. To 
avoid an incorrect combination of cross sections among 
simulated VF/Z+jets samples, only events with the same 
number of reconstructed jets as the number of initial par- 
tons are retained. The background from multijet events, 
in which a jet is misidentified as a lepton, is evaluated 
using the "matrix method" as follows. Two samples of 
VF+jets event candidates are used: a "tight" sample, for 
which the lepton identification criteria are as described 
above, and a "loose" sample in which some of these iden- 
tification criteria are relaxed. The probabilities for true 
leptons to be identified as loose, and jets as tight leptons 
(as determined from independent studies of samples of 
pure leptons and pure jets) yields the fractions of true 
leptons and of misidentified jets in the tight and loose 
samples. 

After the IF-boson and jet selections, we apply the two 



HF-tagging algorithms to the jets. The MC samples are 
normalized to the appropriate luminosity, and corrected 
for differences in HF-tagging and lepton-identification ef- 
ficiencies relative to data. Also, discrepancies in trigger 
efficiency between MC and data are corrected for each set 
of selections. In the following, the e + jets and /i + jets 
samples are combined. Figures ^ and [3] show the exclu- 
sive number of jets in events with at least one SLT-tagged 
jet and at least one SVT-tagged jet, respectively. The 
transverse mass for IF-boson candidate events contain- 
ing at least one SLT or SVT-tagged jet, shown in Fig. 
agrees well with SM expectation. The distribution for 
events with at least one jet tagged with both algorithms 
is shown in Fig. 0] 
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FIG. 1: Exclusive jet multiplicity for VK-boson candidate 
events with at least one SLT-tagged jet. The fourth bin rep- 
resents the sum of events containing four or more jets. 
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FIG. 2: Exclusive jet multiplicity for IF-boson candidate 
events with at least one SVT-tagged jet. The fourth bin rep- 
resents the sum of events containing four or more jets. 

The dominant sources of experimental uncertainty are 
common to both the e+jets and ^+jets selections: (i) a 
6.5% uncertainty on the integrated luminosity, (ii) a 6% 
per jet uncertainty arising from jet-energy-scale correc- 
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FIG. 3: Transverse IF-boson mass for events containing at 
least one SLT- or SVT-tagged jet. 
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FIG. 4: Exclusive jet multiplicity of FF-boson candidate 
events with at least one jet tagged with both the SVT and 
SLT algorithms. The fourth bin represents the sum of events 
containing four or more jets. 



tions and jet identification, (in) a 10% per jet uncertainty 
arising from the HF-tagging algorithms, and (iv) a 10- 
18% uncertainty on the predicted MC cross sections (de- 
pending on sample). The total systematic uncertainty 
on the it background is 16%, 21% on single top-quark 
backgrounds, and 22% on the VF/Z+jets backgrounds. 

No excess is observed in the "doubly-tagged" jet sam- 
ple. We therefore proceed to set a limit on the rate 
anomalous HF-quark production in association with a 
W boson. Because we do not propose a model for such 
production, we do not base this limit on any specific ef- 
ficiency or jet spectrum. We quote limits on the number 
of events beyond SM expectation per exclusive jet bin. 
The 95% confidence level (C.L.) upper limits for addi- 
tional event production in each bin are shown in Table [I] 
These limits are calculated using a modified Frequentist 
(CL a ) method [ll|. 

Assuming that anomalous HF production has the same 
event topology as certain SM process, the above limits 



TABLE I: The numbers of observed and predicted FF-boson 
events with at least one jet tagged by both the SLT and SVT 
algorithms, as a function of exclusive jet multiplicity. Also 
shown are the 95% C.L. limits in the form of additional events. 



Source 


1 jet 


2 jets 


3 jets 


>4 jets 


Data observation 


5 


1 


1 





SM prediction 


5. Oil. 2 


2.0±0.5 


1.0±0.2 


0.4±0.06 


95% C.L. Limit (events) 


6.7 


3.9 


4.1 


3.0 



can be translated into limits on cross sections, 
end, we consider two benchmark scenarios: 



To this 



1. 11 Wbb-like" production in which two b quarks are 
produced in association with a W boson. In this 
scenario, additional light quarks or gluons can be 
produced, and thereby shift the event topology to 
more than two jets. Jets not within the acceptance 
of the detector can also cause the event topology to 
drop to less than two jets. We model this produc- 
tion using efficiencies for SM W/Z + bb production. 

2. "Top-like" production in which a heavy particle is 
produced and decays to a IF boson and a b quark. 
An event can contain two such heavy particles ( Hi- 
like") or one heavy quark ( "single-top- like" ) , with 
additional light or heavy quarks and gluons pos- 
sible for both cases. We model this scenario using 
the cross-section weighted efficiencies for SM it and 
single top-quark production combined. 

We calculate a limit on exclusive jet production for 
each scenario, but first ignore the probability for re- 
constructing the predicted number of jets, providing a 
model-independent comparison of processes with specific 
jet topologies. The remaining efficiency represents the 
effect of VF-boson selection and HF tagging, and lim- 
its for a specific model can be extracted by multiplying 
this value by the efficiency to reconstruct the number of 
jets found in each exclusive jet bin. These results are 
shown in Table [H] To evaluate an upper cross-section 
limit on inclusive jet production for each scenario, we 
reintroduce the efficiency for reconstructing the predicted 
jets. For inclusive PF&5-like anomalous production, we 
sum the first two IF+jets bins, as the contribution from 
the remaining bins is negligible. For top-like anomalous 
production, we sum all VF+jets bins, except the n = 1 
bin, where the contribution is again negligible. Table ITTT1 
shows the 95% C.L. event limits for the combinations 
of jet bins for these two hypotheses, and also the corre- 
sponding anomalous HF production cross-section limits. 
The jet reconstruction efficiency is included in the calcu- 
lations, and the limits contain the expected efficiencies 
for the specified SM processes. 

In summary, we observe no excess beyond the SM pre- 
diction for heavy-flavor quark production in association 
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TABLE II: Cross-section upper limits in pb, based on the hy- 
potheses of "VK&bTike" and "top-like" anomalous production 
of exclusive number of jets. Each value must still be corrected 
for the efficiency of reconstructing the predicted number of 
jets. 

Model 1 jet 2 jets 3 jets >4 jets 
Wbb-like 35.0 9^2 K0 4J5 
Top-like 12.6 8.0 11.3 15.4 



TABLE III: 95% C.L. limits for the number of events summed 
over the indicated jet bins. Also shown are cross-section limits 
based on the hypotheses of "IF&fc-like" and "top-like" anoma- 
lous production for the selected number of jets. 



Source 


1,2 jets 


>2 jets 


Data observation 


6 


2 


SM prediction 


6.9±1.2 


3.3±0.5 


95% C.L. Limit (events) 


6.6 


4.4 


Model 


Wb&-like 


26.4 pb 




Top- like 




14.9 pb 



with W bosons in 164 pb of data in the e + jets chan- 
nel and 145 pb" 1 in the ^+jets channel. Using a sample 
of events containing at least one jet tagged with both 
the SLT and SVT algorithms, we derive 95% C.L. lim- 
its on anomalous heavy-flavor production (Table HJ). Us- 
ing benchmark SM processes, we also derive anomalous 
cross-section limits of 26.4 pb in a Wbb-\ike scenario and 
14.9 pb in a top-like scenario. For comparison, the D0 
collaboration has recently published a similar study in 
the form of a search for Wbb production [J^- Based 
on the two-jet topology, with both jets HF-tagged, that 
study sets a 95% C.L. upper cross-section limit of 6.6 pb 
on Wbb production. 
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